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1.  INTRODUCTION 


The  purpose  of  the  present  work  is  to  re-consider  a well- 
known  resistive  MHD  instability  by  considering  the  effect  of  a 
more  general  form  of  Ohm's  Law.  Specifically,  we  treat  the  m-o, 
finite  wavelength  tearing  mode  in  cylindrical  geometry  for  a 
reversed  field  plasma  equilibrium  t/pical  of  those  which  are 
produced  in  reverse-bias  theta  pinch  experiments?  The  motivation 
for  the  present  study  stems  from  early  experiments  which  showed 
a curious  and  unexplained  absence  of  an  expected  tearing  mode? 
Reverse  bias  theta  pinch  experiments  generally  produce  a closed 
field  configuration  with  reverse  trapped  magnetic  flux;  this 
occurs  due  to  magnetic  line  reconnection  at  the  cold  ends  of  the 
pinch  coil  where  the  resistivity  is  high.  The  formation  of  the 
closed  region  is  also  associated  with  a strong  axial  contraction 
as  the  plasma  attempts  to  find  its  equilibrium.  The  result  is 
an  elongated  (length  to  width  ratio  > 10)  reversed-f ield  plasma 
with  a bounding  closed  separatrix,  which  according  to  the  sim- 
plest classical  picture  should  be  further  susceptable  to  tearing 
modes,  m-o  and  k,  I1  0,  which  would  tend  to  break  the  column 
into  a number  of  cells  each  having  reversed  trapped  flux.  The 
number  and  sire  of  the  cells  as  predicted  by  the  simple  linear 
theory  varies  from  experiment  to  experiment  since  the  linear 
dispersion  relation  depends  sensitively  on  plasma  profiles. 

The  classical  tearing  mode  theory  due  to  Furth,  Killeen  and 

Rosenblut^  carried  out  originally  in  slab  geometry  yields  for 
7 

cylindrical  geometry  the  prediction  that  most  of  the  experimentally 
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produced  reversed  field  theta  pinch  plasmas  should  break  up  into 
reversed  field  cells.  Experimental  results  of  Eberhagen  and 
Grossmann  and  more  recently  by  Linforl  show  that  the  expected 
tearing  mode  does  not  occur  and  the  plasmas  in  question  exist  for 
some  tens  of  linear  growth  times.  Instead  of  the  expected  tear- 
ing a pronounced  rotation  of  the  plasma  column  is  observed;  this 
rotation  increases  with  time  until  the  rotational  frequency 
reaches  a critical  value  at  which  time  the  plasma  goes  catastrophi- 
cally unstable  to  an  m=2  mode.  In  an  attempt  to  explain  the  ab- 
sence of  the  tearing  modes  Kalek^  simulated  the  plasma  behaviour 
with  a linearized  initial  value  MHD  code  assuming  the  plasma 
column  to  be  infinitely  long  with  no  longitudinal  spatial  dependence. 
Kalek's  model  originally  contained  the  resistive  MHD  equations  of 
motion  linearized  about  a rigidly  rotating  plasma  equilibrium; 
results  from  the  numerical  code  showed  that  the  inclusion  of  ro- 
tation predicted  even  larger  growth  rates  than  for  the  static 
equilibrium  case.  The  addition  of  the  Hall  effect  into  the  Ohm's 
Law  yielded  results  which  were  not  inconsistent  with  the  experi- 
mental results,  namely  the  growth  rates  as  calculated  were  at 
least  an  order  of  magnitude  lower  than  what  the  FKR  theory  would 
predict;  this  statement  requires  qualification  in  that  the  rele- 
vant part  of  the  FKR  cylindrical  tearing  mode  dispersion  relation 
is  distorted  in  a manner  to  suppress  the  growth  rate.  The  long 
wavelength  (ka  £ 1)  part  of  the  dispersion  relation  is  affected 
and  a preferred  direction  of  rotation  is  indicated  for  a stabi- 
lising influence.  The  direction  of  the  ion  diamagnetic  drift 
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velocity  is  preferred  for  stability  and  rotation  in  the  opposite 
direction  results  in  enhanced  instability.  Kalek's  numerical 
results  lend  strong  support  to  a beneficial  interplay  between 
the  effects  of  rotation  and  Hall  current  in  the  suppression  of 
tearing  mode  growth  rates.  The  calculation  described  in  the  pre- 
sent paper  attempts  to  analytically  investigate  the  m-o,  k2  r1  0 

tearing  mode  in  reversed  field  plasmas  and  attempts  to  explain 
« 

the  numerical  results  of  Kalek.  Our  results  are  consistent  with 
Kalek's  and  we  find  a stabilizing  influence  of  rotation  combined 
with  Hall  effect;  however,  we  find  that  the  results  are  very  sen- 
sitive to  profiles,  wall  position  and  plasma  parameters  such  as 
density,  temperature  and  rotation  frequency.  The  specific  cal- 
culation presented  in  the  present  paper  deals  with  plasma  profiles 
which  are  representative  of  those  which  are  produced  in  reversed 
field  theta  pinch  experiments. 

It  should  be  mentioned  that  the  use  of  Ohm's  Law  with  a 
Hall  term  has  also  been  used  to  model  the  effect  of  finite  particle 
effects  in  treating  tearing  modes  since  the  Hall  term  is  associ- 
ated with  a coefficient  which  is  proportional  to  the  ion  yy- 

3 

ration  frequency.  In  such  a case  overstable  modes  can  be  found. 

In  the  problem  treated  here  the  m=0  mode  does  not  admit  overstable 
solutions;  this  is  not  explicitly  pointed  out  in  the  text  but  fol- 
lows from  simple  energy  arguments. 

The  paper  is  organized  into  14,  sections.  The  governing 


equations  are  presented  in  sections  2 and  3.  The  governing  equa- 
tions are  non-dimensionalized  in  section  4.  Solutions  of  the 
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governing  equations  in  the  resistive  boundary  layer  are  dis- 
cussed in  sections  5 and  6.  Solutions  outside  the  resistive 
layer  are  discussed  in  section  7 and  matching  across  the 
boundary  layer  is  discussed  in  section  8.  The  tearing  mode 
ordering  with  a small  parameter  is  discussed  in  section  9 and 
10.  Sections  11  through  14  discuss  the  computation  of  growth 
rates  for  specific  equilibria.  A final  section  presents  the 
conclusions  of  the  computations  discussed  here. 
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The  Model 


We  consider  a straight  cylindrical  plasma  column 
whose  equilibrium  profiles  are  considered  to  have 
variation  only  in  the  radial  direction.  The  plasma  is 
considered  incompressible  and  has  a density  profile, 
pQ(r),  that  extends  to  a perfectly  conducting  wall.  The 
equilibrium  magnetic  field,  BQ(r)  has  both  longitudinal 
and  azimuthal  components.  Finally,  the  column  is  sub- 
jected to  a rotation  with  rotational  frequency  fl(r). 

To  study  this  system  we  employ  the  following 
mathematical  model: 


§£  + V* (pv)  = 0 


dV 

p + pv • Vv  = -Vp  + JXB 


SB 

^ + 7xE  = 0 


J - — VxB 


-VXB  + nJ  + <J*B 


7 • v = 0 


7*  § = 0 
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where  v is  the  plasma  velocity,  § the  magnetic 
field,  p the  plasma  pressure,  J the  current  density, 

£ the  electric  field  and  n»  e are  resistivity  and  Hall 
current  parameters  respectively  which  are  taken  to  be 
constant.  The  mathematical  model  specified  by  Equations 
(1)  - (7)  differs  from  Ideal  MHD  only  through  Ohm's 
Law  which  includes  the  effect  of  finite  resistivity  and 
the  Hall  effect.  Our  aim  in  this  report  is  to  study 
the  time  behavior  of  small  perturbations  about  an 
equilibrium  specified  by  vQ(r),  BQ(r),  pQ(r)  and  pQ(r). 
The  perturbed  velocity  field,  magnetic  field,  density, 
pressure,  current  density  and  electric  field  are  repre- 
sented by  u,  b,  Sp,  5p,  j and  5E  respectively. 

The  linearized  version  of  equations  of  (1)  - (7) 

are: 

<Sp  + V*  (5pyQ)  + V*(pQu)  = 0 (la) 

It  PQ  « ■ “*P*  + 2*?o  + JQxb  - PQ  yQ  (2a) 


- p0  yo-vy  “6py0*vY0  + J-?(§*b) 


d . 

yt  * = 

-Vx6E 

(3a) 

1 - J-  7x* 

(4a) 

uE  — -uxo  - V xu  + i| j 

+ i — { j X D + vj  xjj  - 

- — o o ) v Da ) 

- - ~o  -o  - i 

p J -o  -o  - 

p -o  -o 

" ■ "" 
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V*v?  • 0 (6a) 

- 0 (7a) 

The  starting  set  of  equations  for  the  computation 
can  be  obtained  from  system  (la)  - (7a)  by  Fourier 
transforming  the  perturbed  quantities  in  the  ignorable 
Q,  s directions  and  time,  t.  As  mentioned,  we  consider 
only  the  mode  m « 0,  k finite  where  m is  the  asimuthal 
wave  number  and  k is  the  axial  wave  number.  Thus, 

• &(r)  exp(i[ut  - kz]  ) 
u « u (r)  exp(i[uit  - kz] ) 

etc. 


For  ease  of  notation  we  will  omit  superscripts  from 
the  perturbed  transformed  variables. 

After  transforming  equations  (la)  - (7a)  and  rear- 
ranging we  arrive  at  the  following  sixth  order  system  of 
ordinary  differential  equations: 


,’i  «(.rV 

o r 


3r  (pow  r IF 


, . d (rb  ) , , dp 

1 ~ cf?  (f  ■^rr~)  k po(n>  337 
0 po 


-w2  ♦ 4>12  ♦ 2r>7  ^Kr  + i(^(-  J g - *Vl 


ihhj  )b  . * !!i  k2 

o,  a r 1 wp0  r 


> b. 


to) 


(8) 
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3.  The  Equilibrium 

We  consider  an  equilibrium  magnetic  field  profile 
with  vanishing  Bq  z on  a cylindrical  surface,  r - rQ,  in  the 
plasma,  i.e.#a  reversed  field,  while  the  component,  3O,0 


of  the  field  is  taken  to  be  finite  on  r - r^. 

o 

Since  the  plasma  is  taken  to  be  incompressible,  we 
may  impose  an  arbitrary  density  profile.  Finally,  the 
equilibrium  flow  is  taken  to  be  a rotation  with  variation 
in  the  radial  direction,  and  no  flow  in  the  axial  direction. 

The  above  equilibrium  may  be  summarized  as: 


So‘r>  • Bo,e  * * Bo,z  z ' 


b„  . <*„>  ■ 0 
O f z o 


vQ(r)  - rfi(r)  0. 


PQ(r)  arbitrary. 


4.  Dimens ionalizat ion: 


We  dimensionalize  with  respect  to  the  equilibrium 

values  at  the  neutral  surface  r - r^  as  follows: 

o 


2 

oj  a u 
o o 


c df  . . 

Po^O  dr  ° 
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P0(r) 

,Q " >7*7 


f f df 

' * adf  (r . ) ' dr  “ < TF, 

a?  ° a?( 


df 


dr 


<rol 


d2f 

7?^ 
a 


ar(ro> 


r-r. 


/ T 


dp 

D - - 2 -2  , B.  (r  ) 

s dr  u o 

o 


E B 


0/6 


a2  df.  , 


E J. 


J_  « . Y “ “ 


Z _df(l  ' 
o 


aa?(r-> 


fi  - — , k = ka  , ~ * — 

u>  dr  u> 

o o 


u> 


(polJo) 


a df  lbi 

“IT?  a?  (ro)  ' qr  '_rdT 


a d?(ro} 


ikb. 


qe  " _2“d7 


a a7tro> 


With  these  dimensionalizations : 


T(t) 


and 


►t  a a 

aia  * o n m 


-li- 


ne rewrite  equations  (8)  - (10)  in  dimensionless 
form  for  a layer  centered  about  rQ  that  is  sufficiently 
small  to  neglect  variations  in  the  equilibrium  quantities. 
This  layer  will  be  referred  to  as  the  boundary  layer, 
and  the  equations  governing  the  behavior  of  the  resistive 
layer  will  be  known  as  the  resistive  boundary  layer  equations 
or  RLE.  To  leading  order  in  x and  omitting  the  super- 
scripts on  the  dimensionless  quantities,  the  RLE  are: 


-(i  - IcVc^  - -y-)qr  ♦ (Si  - -^)q. 


A,. 

d2q 

- r r 

dx2 

T 2 
dr^ 

X 

i 

^|M0 

1 

d2qr  _ 
dx2 

-YT£r  + i 

d^0 

d2q 

r 

dx2 

ex  2 

dx^ 

->  r D 

i‘n  O S , 


■ondfi,  , ,2  2,  _ 


. ,nx  2iBe, 


= “YT£r  + exq0  + (y+k  n + iea2x)qr 


■2  * r * icvi " Dsnq9 


r jp  2iB_ 

+ teDsTU+a2Be]  + fix  [1  + + Y(a2x  - ^ 

o 


.2,  10t2  2 , ro  d«, 


(eta,  - + Ds  + k‘x]  - -f  t-  + ^^)qr 


(10a) 


Analysis  of  the  RLE 


Scaling: 

We  wish  to  study  the  interesting  behavior  governed 

by  the  RLE  in  the  limit  of  small  resistivity  and  Hall 

parameter  n and  e.  Following  the  suggestion  of  reference 

(9],  we  adopt  the  following  scaling  based  on  a resistive 

layer  thickness  of  the  order  of  the  skin  depth  with  respect 

2 

to  the  diffusion  time  scale,  i.e.,  yx  - n where  x y. 

y * 4>s  , x = , n = 4>3n  , = 0(1) 


e * 4>e 

**  5 + • • • » = tpty  + • . • , 

qQ  = ? + *.. 


Also  denote: 


P 


fo  dfi 
dr  ' 


r 


o 


where  the  superscripted  variables  are  considered  to  be 
0(1)  quantities  while  $ is  a small  scale  factor.  The  scaling 
of  the  Hall  parameter  and  angular  velocity  were  chosen  so 
that  terms  with  e and  ft  would  neither  dominate  nor  disappear 
from  the  scaled  equations  in  the  lowest  order. 


Finally,  we  again  omit  the  superscripts  in  writing  the 
scaled  equations: 


i2  ♦ x ■ n2(l  - J2  Ds  + PK  + (—  - i6)c 


H — " ~S*£  + exy  + Sip 
dx2 


- - (S  + + iefi  [ 1 - 

dx2  b 2 


♦ (xflfl  + P]  - i6S)i  - P(hp 


Making  use  of  (12),  equations  (11)  and  (13)  can  be 
rewritten: 


dx2 


(ifs  + fl2{i  _ fo^S  + pJ)  ^ (£x  _ ex; 

n * s n 


- id)  C - <P 


(11a  ) 


n -V 

dx2 


x2  r2  2 * 

(S  ( 1 ♦ §-]  ♦ {j-  x2  + ic3  (1  - -§  Ds])c 


♦ (xfl ll+P]  - ifiS  - S)£  + (^  - pfl)ij.  . (13a) 


Equations  (11)  - (13)  show  that  in  the  presence  of 
azimuthal  components,  Bg,  of  the  equilibrium  magnetic  field, 
overstable  modes  are  possible  due  to  the  imaginary  terms 


in  the  equations.  For  this  analysis  we  suppress  these 
terms  by  considering  only  long  wavelength  modes  with 
sufficiently  small  k which  results  in  the  elimination 
of  6 and  a2  from  equations  (11)  - (13)  to  leading 


order 
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6. 


Asymptotic  Behavior  of  RLE: 

The  system  (11a),  (12),  (13a)  constitutes  a sixth 
order  homogeneous  linear  differential  equation  with  an 
irregular  singularity  at  x » ®.  The  system  may  be 
rewritten  in  a standard  form  for  irregular  perturbation 
analysis  (see  reference  (10)): 


du 

- x A (x)  u 

which  has  a set  of  solutions  of  the  form: 

X x2  D 

« * • • • «,  A 

u “ c e X z ( x ) 

where  c is  a constant  matrix  and  X , D are  diagonal 
matrices. 

The  elements  of  X determine  the  exponential  nature 
of  the  asymptotic  solution  to  the  inner  equations.  These 
elements,  X^,  are  the  eigenvalues  of  the  matrix, 


A («*>) 


0 0 0 1 
0 0 0 0 
0 0 0 0 


Jo  - ^ 0 0 


eQ 

p. 


0 


0 


0 0 0 0 


0 0 
1 0 
0 1 

0 0 

0 0 
0 0 


Thus, 


1*2 


0,0 


3,4 


(ns) 


1/2  2 


2 1/2 
d * (1  + r)  ) 


1/2 


5,6 


(nS) 


rm  <1 1 (1  + f->  » 


where  r4 


c2s 


At  this  point  it  is  not  clear  whether  A^» 
k ■ 3, 4, 5, 6 have  pure  or  oscillating  exponential  growth 
or  decay,  or  whether  they  have  pure  oscillatory  be- 
havior as  x + ».  In  fact,  we  can  show  that  only  if. 


S 


0- 

2 2 
(ot  +1)  e 


(-1 


for  o an  arbitrary  real  number,  can  A^  be  pure 
imaginary  in  which  case 


A 


k 


Thus,  all  these  solutions  either  decay  or  grow 
exponentially  as  x + » in  the  class  of  solutions  with 
positive  growth  rates,  i.e.,  Re  s > 0 

Finally  we  consider  the  two  algebraic  equations 
corresponding  to  A1  and  \^.  To  determine  the  asymptotic 
behavior  of  these  solutions  as  x we  consider  the 
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dominant  terms  from  equations  (11a)  and  (13a),  as 
x •: 


. £i  ; sj 

n n xn 


(-  - — ) £ + (i  + )y  - -SH  y 
x n s n T xn 


Solving  to  leading  order  in  x for  £ and  y yields: 


£ = * , 
* x 


Y 


sn  . 
T * 

x 


(14) 


(15) 


Thus, 


Y 


x * 


(16) 


Also  from  (12), 


djt  s d2(x£) 
dx2  dx2 


- « s 

dx^  b 


(17) 


Replacing  (14)  and  (15)  in  (11)  and  noting  that 


x df2  (xC)  = 3x  (*2  if 


-19- 


7.  Solutions  Outside  the  Resistive  Layer; 

The  behavior  of  the  plasma  outside  a narrow  region 
surrounding  the  neutral  surface  is  governed  to  leading 
order  by  ideal  MHD,  i.e.,  c and  n can  be  taken  to  be  zero. 

In  terms  of  the  same  dimensionalizations  as  the  resistive 
layer  equations,  the  outer  equations  reduce,  near  marginal 
stability,  i.e.  Y 2 0,  to  the  following  second  order  differ- 
ential equation: 


d ,f2  d 

3F(F"  3F 


<r£r))-f2k2£r 


po 

S 


rr‘2(V§  aX 


(19) 


Near  rQ  equations  (11)  - (13)  reduce  asymptotically  to: 


d 

a? 


2 d*r 

a ^ 


o r«Dc 

-fi2(_oS  . P)r 


(20) 


which  has  solutions. 


€ 


T 


1-h 


(i  - a2(roDs  - p))1/2. 


(21) 


where  h 


(21a) 
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For  the  purposes  of  numerical  integration,  it  is  con- 

er 

venient  to  set  W - which  transforms  (19)  to, 

^r 

3F" 


♦ P(x)  W2  + q (x ) W - 1 


where. 


W(0)  - W(l)  - 0 


poo  . i * m*  ♦ k*  . V£l 

x*  4f ^ 


and. 


q (x) 


- (2  df/dx  + 1 
f x' 


From  (22)  and  the  asymptotic  expressions  fix,  = 1» 
we  find  that  the  criteria  for  entering  the  asymptotic 
regime  governed  by  (20)  is  that 


f1  ,2  I « 1 • 

o s 


In  this  formula,  the  value  of  x such  that  M-.  01  may  he 
taken  as  a measure  of  the  thickness  of  the  boundary  layer,  i.e., 
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8. 


! 

! 


b.£. 


.01 


x2dsq2 
o s 

8 


(23) 


We  note  that  according  to  (23)  the  boundary  layer 
thickens  as  the  plasma  rotates  faster. 


Matching  Inner  and  Outer  Solutions: 

In  order  to  match  the  outer  solutions,  the  two 
exponentially  growing  solutions  to  the  inner  equations  must 
be  set  equal  to  zero,  leaving  the  two  algebraic  solutions 
for  matching. 

The  condition  that  the  algebraic  inner  solutions  match 
the  outer  ideal  MHD  solutions  can  be  expressed  as: 


outer  inner 

lim  Cr(t)  = lim  £ (x)  (24) 

T+0  x->» 

From  (18)  and  (19)  it  is  clear  that  the  inner  solutions 
can  be  matched  asymptotically  to  the  outer  solutions  in  the 
sense  of  (24) . 

We  will  be  investigating  a class  of  resistive  eigenmodes 
called  tearing  modes,  that  are  sensitive  to  the  matching 
conditions  between  the  inner  and  outer  regions.  We  subdivide 
the  plasma  radius  into  three  regions.  Region  I is  the  outer 
region  closest  to  the  magnetic  axis.  Generally,  when  treat- 
ing boundary  layers,  the  boundary  layer  is  called  the  inner 


ever,  we  have  two  outer  regions,  one  of  which  is  in  some 


sense  inside  the  boundary  layer.  We  can  just  call  it  region 
I.  Region  II  is  then  the  inner,  resistive  boundary  layer, 
and  region  III  the  other  outer  region. 

From  Equation  (21)  we  have  seen  that  the  solution  in 
Region  I near  the  boundary  layer  approaches  the  form: 

5r,I  = AI  lT|h  + Bi I T | ~1-h  (25) 

Likewise  in  Region  III, 


«r,III  = AIIlMh  * BIIlltr1‘h 


(26) 


Similarly  in  Region  II, 


5 = 


An,i>ih  + as  x * - 


AII,rI xl  ” + BIT.plxl'1'h  as  x * * 


II,  R1 


The  matching  condition,  (24)  can  be  expressed  as, 


B. 


It  2h+1  AIIiL 


(?> 


B 


(27) 


II, L 


and. 


‘III 


B 


III 


.1,  2h+1  AI I , R 

bii,r 


(28) 
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2 1/2 

for  2h+l  ■ (1  - £1  (r  D_  - P) ) ' >0  and  where  $ is  the 

o S 

small  scale  factor. 

A A 

We  expect  for  most  cases  that  _I  and  III  will  be 

B BIII 

finite  numbers.  Thus  to  obtain  matching,  the  growth 
rate  must  be  adjusted  in  the  resistive  layer  so  that 
Ajj  L - 0.  The  exceptional  case  occurs  when  Bj  vanishes. 
According  to  Newcomb's  stability  criterion  [11],  this 
is  the  condition  for  marginal  ideal  stability.  It  is  not 
surprising  that  a system  that  is  close  to  being  unstable 
in  the  ideal  theory  will  have  a strongly  modified  resistive 
instability. 

These  matching  constants  can  be  calculated  numerically 
from  the  outer  layer  equations  as  the  limit 
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In  terms  of  W,  (29) can  be  rewritten. 


III  _ lim  - 


Tf2h_1  h+1  + J 

1 T 1 

-h  ♦ * 


(30a) 


5*-  - lim  - 
BI  x-»-0 


i -2h-l 


The  quantity  matched  to  the  inner  layer  is, 


(30b) 


Our  objective  in  this  report  is  to  gain  some  infor- 
mation about  the  eigenvalue  problem  by  purely  analytic 
means.  To  do  this,  we  consider  the  growth  rate, 

S,  and  h - 0.  It  will  prove  convenient  to  match  the 
r-component  of  the  perturbed  magnetic  field  at  the  edge 
of  the  resistive  layer  rather  than  the  plasma  displacements, 
in  which  case  the  asymptotic  solutions  are 


q„  T “ Br  + ATt  • as  t -*•  0 

i f JL  1 J. 


<lr.HI  = Bttt  + Atttt  ' as  t 0 • 
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and  , 


* " 


IIf  L + AII,L 

bii,r  + aii,r 


as  x-*--00 


as  x 00 


As  in  (27)  and  (28)  the  matching  conditions  become. 


<P  B 


II, L 
II, L 


0(1) 

BII, 


(30) 


We  have  seen  that  the  behavior  of  the  exponential 
solutions  of  the  RLE  are 


5 


exp 


{± 


kx 

2(nS) 


172} 


for  k constant  , 


as  x Thus  as  S becomes  small,  we  must  consider  the 

1/4 

length  scale  to  become  small  as  S , otherwise  the 

exponential  solutions  would  grow  infinitely  rapidly  making 

the  equations  badly  posed.  That  is,  there  is  an  inner 

layer  of  significant  behavior  that  is  even  thinner  than  the 

1/4 

resistive  skin  depth,  by  another  factor  of  S . 

We  find  the  interesting  behavior  within  this  new  scaling 
by  expanding  the  parameters  and  variables  of  equations 
(11)  - (13)  in  powers  of  S as  follows: 


€ * eo  + ••• 


<l>  * S1'4  <PQ  + s7/4  ^ + • • • 


c * s 


3/2  y + ... 
1 


x - sl/4  x , e - -J72  e ' = S3//4  fi 


Equations  (11)  - (13)  reduce  in  the  zero-th  order  to 


d2i 1>C 
dx2 


= 0 


<*2 


x2  —2  roBS  £ ~2  v 

<7  ' a ‘“H  + (0x  - t x ' <0*^0 


«2t„  -2  l2  — lx2  - x 

— 5s  " x I1+f  l«o  + + (c  f - P0)»o 


dx 


In  the  first  order: 


a2*, 

" 7=2  * -x«o  + £x;o  + *o 


(0)  (0)_ 

^o  = BII  + AII  x 


Ola) 


(31) 


Thus 


But  since  A - 0 for  matching,  we  are  justified  in 

(0)  11 

setting  AXI  ■ 0 in  the  zeroth  order  and  accounting  for 
the  x dependence  in  the  next  order.  In  fact,  from  (31), 
and  in  view  of  (14)  and  (15),  we  have 


d2i(» 

£ - 0 as  x -*•  ±“ 

dir 


which  has  solutions. 


*1 


(1)  _ 

aii,l  x as  x"-°° 


(1)  - 

aii,r  x as  x - - 


where  the  constant  solution  has  been  incorporated  in  the 
zeroth  order  solution.  Thus  the  total  solution  is. 


- 7/4  1/4  (0) 

<(>(x)  - s'  AXI  x + S 7 BII; 


And  since. 


di{/ 

dx 


(x) 


1 difr 
S1/4  dx 


(x) 


making  use  of  (25),  (26)  and  (31),  we  arrive  at  the  dis- 
persion relation, 


" 


28- 


(33) 


(34) 


(35) 


9.  The  Tearing  Mode  Equations: 


From  the  original  scalings  of  section  5,  and  equation 
(34),  the  tearing  mode  scalings  can  be  rewritten  in  terms 
of  the  scale  factor  6 as 


n a 5*  □ 3*  dO  q d A 

n = en,y  = 0Q,5?=e3F 


x = 02 *x  , e = 0e  , 0 = 0A 


qr  = q2^0  + 04^1  + 


5 = 5 + ... 

r 


(36) 


qe  = 9 ?o  + * ‘ * 


P 


and  letting  P * ^ gp 


eDS 

■\nd  q - — . 


(36a) 


The  following  tearing  mode  equations  result  from  scaling 

the  basic  equations  (11)  - (13)  according  to  (36), 

r D 

omitting  superscripts  and  considering  |P  - ° - | <<  1 


2 

m (£_  [l  - -2-S  + P]  + jMc  + (A|  - 

Q2  4 Qn  o q3  Q2n  o 


dx‘ 


a 1 
1"  r 


+ (- 


nQ 


)6 


(37) 
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2 2 


— 7 - <r  I1  + ' ^2QH o + (iQ  + ^ 


+ <’  >1^  -Vql  +^,e 


(38  ) 


d2^ 


-li  - - 22  5 + «c  + 2 g 
dx2  n n '•o  n 


(39  ) 


where  the  superscripts  have  been  omitted  from  the  scaled 
variables. 

We  study  the  asymptotics  of  system  (37)  - (39)  by 
expanding  the  independent  variables  as  follows: 

(40) 

(41) 


Replacing  these  expansions  in  (37)  - (39) , we  obtain 
the  result  that 


6 , 


C2  " 


-6QA  , e. 


[g  (i-  - c^)  - eA] 


(42) 


and* 

• 2r  D 

- 8(1  + - — In  x 


...]  + [c2  + 8(01^-  — ) In  x)x 

o 


where  0 and  are  chosen  to  match  the  outer  solution • 
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In  order  to  study  system  (37)  - (39)  analytically 
we  are  motivated  by  equations  (42)  to  make  the  following 
change  of  variables: 


With  these  substitutions  (37)  - (39)  can  be 

r D<; 

rewritten,  assuming  that  |— —■  - - P|  <<  1: 

2 F D 

— 2 * [A2  (1  - — ° -)  + z2K  + A [ Az  - rz2]c  + (a^ 

dz 


z 


(43) 


(44) 


t 

' 





• m I 
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dz2 


- (zA  (1  + q)  - rz2) 5 + z2 (1  + r2)l 


* ("  - 

A Q"’A 


3747  (“l  ' ’ 


(45) 


"4  * T7I  (1  ' + r^> 

dz  n 


(46) 


In  the  text  that  follows  we  will  simplify  the 

computation  by  taking  a.  = ~ . . 

x o 

The  dispersion  relation  (32)  - (34)  can  also  be  rewritten 
in  terms  of  these  substitutions: 


1 d£ 
6 3x 


jV4 

.3/4 


where , 


. f5 

un 
♦*00  ' 


lim 

Z+oo 


(1  - z£  + rA?)dz 


-z 


and. 


1 dto 
6 dx 


1 lim  rd^  . , d£  . n 
6 X—  dx  (X)  dx  ' xn 


Thus,  the  dispersion  relation  is: 


(47) 


(48) 


. il 
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and  making  use  of  (36)  we  have: 

9 - n3/s  4')  4/5 

where  the  dimensionless  quantities  y and  n have  been  used 

in  place  of  the  scaled  variables  Q and  n , and  P depends 

on  the  parameters  in  equations  (44)  - (46)  » namely, 

q,r,A,  and  D . 

s 

A Special  Case: 

Our  objective  in  this  section  is  to  find  a special 
class  of  analytic  solutions  to  the  inner  layer  equations 
(44)  and  (45),  and  use  these  solutions  to  evaluate  the 
dispersion  relation  (48)  explicitly. 

The  case  we  consider  is  characterized  by, 

r D 

-2-2-  « l,  r <<  1,  P = 0,  q - - 1,  A * 0(1) 

This  corresponds  to  a Theta  pinch  equilibrium  with 


(49) 


rigid  body  rotation  and  density  gradient  opposite  in 
sign  to  the  angular  velocity. 
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V. 

From  (44)  we  see  that  only  the  odd  part  of  £ 
and  C contribute  to  the  dispersion  relation.  Furthermore, 
in  the  limit  of  small  Q,  the  second  term  in  (44)  can 
be  treated  as  a perturbation,  thus  only  the  odd  part 
of  X must  be  determined  in  order  to  derive  the  dispersion 
relation.  But,  referring  to  (43)  and  (44),  and  using 
the  result  of  the  asymptotics  summarized  in  (14)  and  (15), 
all  terms  involving  r may  be  treated  as  pertubations  and 
neglected  to  leading  order.  Thus  (44)  - (46)  reduces  to: 


(50) 

1 

<iz 

(51) 

1 

& 

1 r o3/2 

T74  Vi  (1  - za  dz 

(nQ)  ' nx/ 

(52) 

where  we 

denote  the  value  of  A corresponding  to  the 

special 

case 

q ■ -1  by  We  let. 

i - t**  /J  I b2k  H2k(*) 

(53) 

« ■ t <*>  • 

(54) 

-35- 


and  make  use  of  the  identities, 


1 - /2  £ 


z - /2  c 


-*72 


H (z) 
2k 

22kr (k+1) 


-z2/2  H2k+l(z) 


22kr (k+1) 


Solving  for  £ we  find. 


12  + 


Ao(4k+l 


4k+5 


)] 


2k+l 


2ZK  2T (1+1) 


(4k+3  + A2 
o 


Replacing  £ in  (51) , and  using  the  orthogonality 
relations  for  Hermite  polynomials  yields, 

„ ;VI  . 

4 ' PtV 


where, 

3+A2  . 

-2  r(~l^  -2  T{V 

F(A  ) - -4(1  + A‘)  ry~  + 2A‘  — 4~ 

° O „ Or, 3, 


5+a: 


r(- 


:) 


r(4) 


Growth  rates  can  be  obtained  by  solving  the  dispersion 
relation,  A ' * A,  for  the  growth  rate  as  given  by  (49) 
where  the  quantity  A'  is  determined  by  integrating  the  ideal 
equations  of  motion  forward  from  r * 0 and  backward  from 
r " rw(rw  = wall  position)  and  then  forming  the  logarithmic 
derivative  as  indicated  in  Eq.  (30a).  From  (49)  we  notice 


that  the  stability  criterion  is  i'  > 0 , that  is,  for 
positive  values  of  A*  the  plasma  is  unstable  and  positive 
growth  rates  can  be  calculated.  For  a'  0,  which  depends 
only  on  equilibrium  quantities  outside  the  resistive 
layer,  the  plasma  is  stable.  We  note  that,  in  general, 
the  dispersion  relation  will  have  to  be  solved  iteratively 
since  from  (43)  » AQ  depends  on  the  growth  rate. 

11.  Validation  of  the  Analysis 

The  ultimate  goal  in  this  work  is  to  compute  the 

stability  (growth  rates  and  stability  boundaries)  of 

realistic  linear  plasma  configurations  to  the  tearing 

modes  considered  here.  Before  proceeding,  it  is  impro- 

tant  to  validate  the  analysis  with  a specific  example. 

Let  us  consider  the  special  case  of  the  last  section. 

We  will  show  by  a numerical  integration  of  the  governing 

equations  in  the  resistive  layer  that  i)  the  asymptotics 

are  correct  and  ii)  the  analytic  solutions  given  by  (53) 

and  (54)  are  correct.  We  use  a technique  for  solving  Eqs. 

(50)  “ (52)  which  converts  the  boundary  value  problem 

into  an  initial  value  problem  [12].  This  technique  has  been 

used  successfully  elsewhere  in  a number  of  resistive 

boundary  layer  calculations  [13].  We  have  chosen  the  special 
2 

case  q * -1,  A = 0.75  and  have  numerically  integrated  the 
resistive  layer  equations  for  £ and  y as  a function  of  z 
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for  0 <_  z <_  4.  Shown  in  figure  2 is  a curve  of  £ vs.  z; 

the  crosses  denote  the  solution  as  given  by  the  analytic 

expression  (89) , and  the  circles  denote  the  results 

from  the  numerical  integration  of  Equation  (50)  and  (51). 

Figure  2 presents  a similar  comparison  between  analytic 

and  numerical  solutions  for  y vs.  z.  The  asymptotic 

expansions  for  | and  y for  large  z,  i.e.  t.  - 1/z  and 
_ 2 

y - 1/z  , are  borne  out  by  both  the  numerical  and 
analytical  solutions.  Further,  comparison  of  the  solu- 
tions for  the  entire  z range  indicates  excellent  agree- 
ment. We  therefore  believe  that  the  numerical  technique 
for  solving  the  RLE  equations  is  satisfactory. 


r 


12.  Effect  of  Hall  current  on  the  growth  rate: 

From  equation  (49)  we  see  that  the  effect  of  Hall 
current  on  the  growth  rate  is  entirely  contained  within 
the  quantity  F(A,q,r,Ds).  The  value  of  F is  governed, 
for  the  most  part,  by  the  parameters  of  the  inner  layer, 
namely,  e,o  and  the  local  values  of  the  angular  velocity, 
Q,  and  the  density  gradient,  Ds.  Thus  the  quantity  F 
maybe  considered  as  a criteria  of  merit  which  gives  us 
a means  of  evaluating  the  effect  of  Hall  current  on  the 
growth  rate  of  the  tearing  mode  as  compared  to  growth 
rates  in  the  absence  of  Hall  current.  For  example  in 
the  absence  of  Hall  current  and  the  case  of  rigid  body 
rotation,  e * q * 0 and  F = F(A,0,0,D  ).  Thus  larger 

S 

values  of  F in  the  presence  of  Hall  current  compared  to 
no  Hall  current  would  according  to  (49)  tend  to  depress 
the  value  of  the  growth  rate  whereas  smaller  values  of  F 
would  tend  to  have  the  reverse  effect  of  enhancing  the 
growth  of  the  tearing  mode. 

It  is  important,  therefore,  to  examine  the  criteria 

of  merit,  F,  more  closely.  We  will  demonstrate  later  that 

for  realistic  values  of  the  parameters,  F depends  only  on 
-2  - —2  - 

A and  q.  In  other  words  F “F(A  ,q)  , and  the  parameter  r 


along  with  the  explicit  dependence  of  F on  D through  the 
r D s 

term  in  equation  (44)  may  be  ronsidered  as  small  pertur- 


bu ticns  • Tbi  3 3 1 ir.p i 1 1 i c i 
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siderably  and  enables  us  to  illustrate  the  dependence  of 
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F on  A and  q by  Figure  3 . If  we  consider  rigid  body 

rotations,  the  curve  q = 0 corresponds  to  the  case  of  no 

Hall  current.  From  this  figure  it  is  evident  that  values  of 

q > 0 result  in  enhanced  stability  to  tearing  modes  whereas 

values  of  q < 0 result  greater  instability.  For  the  case  of 

- eDs 

rigid  body  rotation,  from  ( 43) , q = . Thus  rotation  of 

the  plasma  in  one  direction  relative  to  the  value  of  the 
density  gradient  at  the  neutral  point,  enhances  the  stability 
to  tearing  modes  while  rotation  in  the  opposite  direction 
results  in  greater  instability.  This  result  was  observed  by 
the  computer  analysis  of  Kalek  [7]. 

Figure  3 also  shows  that  the  effect  of  rotation  on  F, 
and  therefore  the  growth  rates,  is  comparitively  small  for 
the  case  of  no  Hall  current  due  to  the  relative  constancy 
of  the  curve  q = 0 = e compared  to  curves  with  larger  values 
of  q. 

It  should  also  be  noted  that  in  the  case  of  no  Hall 

r - 

current  but  a sheared  velocity  profile,  q = — ° ^ , and  once 

20  dr 

again  Figure  3 shows  that  sheared  velocity  can  also  effect 
the  stability  of  tearing  modes.  However,  this  point  will 
not  be  pursued  in  this  paper.  Hereafter,  we  will  special- 
ize our  analysis  to  rigid  body  rotations. 
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13.  Effect  of  the  equilibrium  and  position  of  the  outer  vail 
on  the  tearing  mode: 

The  equilibrium  determines  the  parameters  of  the  inner 
layer  and  thus  influences  the  tearing  modes  through  the 
criteria  of  merit,  F.  The  equilibrium  also  affects  the 
tearing  modes  through  the  outer  equations  via  the  parameter 
A'.  Furthermore,  since  by  equation  (49),  the  value  of  a’ 
is  extremely  sensitive  to  the  position  of  the  outer  wall 
sur rounding  the  vacuum  region,  we  see  that  growth  rates  and 

_ 9 

thus  values  of  the  inner  parameter,  A defined  by  (43),  are 
very  sensitive  to  the  outer  wall  position.  Thus  values  of 
A determine  not  only  the  effect  of  rotation  but  also  of  the 
position  of  the  wall  on  the  value  of  the  criteria  of  merit, 
F. 

We  now  show  how  the  Hall  current  rigid  body  rotation, 
and  wall  position  effect  the  growth  rates  of  tearing  modes 
for  a family  of  rigid  rotor.  Theta  pinch  equilibria.  To  be 
exact,  we  will  determine  an  overestimate  of  the  ratio  of 
growth  rate  with  Hall  current  to  growth  rate  with  no  Hall 
current,  and  show  how  these  ratios  are  affected  by  angular 
velocity  for  fixed  positions  of  the  outer  wall.  The  next 
section  is  devoted  to  describing  the  equilibria. 


14.  A family  of  equilibria: 

In  this  section  we  apply  the  theory  to  a class  of 
reversed  field,  isothermal,  rigid-rotor  Theta  pinch  MHD 
equilibria.  We  assume  that  all  the  mass  and  rotational 
energy  i»  carried  by  the  ions  and  that  the  family  of 
equilibria  is  parametized  by  the  angular  velocity,  ft, 

N 

and  characterized  by  a constant  number  of  particles,  L. 

The  equilibria  are  one  dimensional,  and  although  their 
domain  of  definition  is  from  the  plasma  center  line,  r = 0, 
to  r = ®,  the  density  is  effectively  zero  for  r >_  a;  where 
a is  taken  to  be  the  plasma  radius.  Furthermore  all  equi- 
librium quantities  are  constant  from  the  plasma  edge  to 
the  wall  of  a vacuum  region,  rw,  surrounding  the  plasma, 
and  the  plasma  column  has  length  L and  is  topologically 
homologous  to  a torus.  Figure  1 illustrates  this  configuration. 
The  equilibrium  is  characterized  by  the  force  balance. 

If  (p  + ij!r  b2)  = n2rp  <55) 

and  the  equation  of  state, 

P - cp  , c « ^ (56) 

mi 

where  K is  Boltzmann's  constant,  itv  is  the  mass  of  an  ion,  ar.d 
T is  the  plasma  ion  temperature. 


■J 
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The  number  of  particles  (ions),  n , is  determined  by, 

L 


nl  - f 
"i  J< 


P (r) rdr 


and  expressed  in  terms  of  a characteristic  length,  r^,  as 


N it  , , 2 
L 0,o)ri 


We  consider  the  following  ansatz. 


p(x)  ■ p(o)  sech  ( 


- «o> 


and, 


b2  - B«o  (3  + d-'s)  tanh(  J-]2  - £o)j 


where,  all  lengths  are  normalized  to  the  plasma  radius, 
i.e. , x - ^ . 

A 

Replacing  (58)  and  (59)  in  (55)  and  (56)  and  character- 
izing the  ratio  of  the  thermal  to  the  magnetic  energy  by, 

„ p(o)u 

$ - 2c  o ' 


and  the  ratio  of  the  rotational  to  the  thermal  energy  by 


t 
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yields, 


n*2  _ 26  (1-6)  sech  KQ 


P ■ (1-6) 2 sech2£ 


x2  . x2  sech  {D 
p 1+tanh  E 


The  equilibrium  is  therefore  entirely  specified  once 
the  physically  meaningful  variables,  B^,  p (o) , T,  NL  and  ft 
are  specified.  However,  for  our  purposes  it  is  more  con- 
venient to  specify,  Ba,  N^,  5,  and  x^. 

In  terms  of  the  latter  set  of  parameters,  the  plasma 
temperature,  rigid  body  rotion  and  the  neutral  point,  xq , 
of  the  plasma  are  given  by. 


- 2 7i 6 (1-6)  (1+tanh  £ ) 

l*  * — 

m.  u Nt 
1 O L 


where. 


c - fl2a2x2  (1-tanh  ^)  i±f-) 


1 
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and, 

xQ  - xl  (1-tanh  50)1/2  (C0+|'ln(l-2y))1/2 

The  object  of  this  study  is  to  show  the  effect  of  the 
Hall  current  dominated  boundary  layer  on  the  growth  rates, 

Y,  of  the  tearing  modes  expressed  by  the  dispersion  relation, 

5 . *3/5  Ul^_)4/5  <«> 

where,  once  again,  the  superscripted  variables  are  normalized 

with  respect  to  the  boundary  layer  quantities  and  the  variables 
-2 

q and  A governing  the  inner  equations  are  defined  by  (36a) , 

and  (43).  Values  of  the  mismatch,  a',  are  obtained  by 

integrating  the  outer  equations  as  described  in  section  7 
-2 

while  F (A ,q)  is  obtained  by  integrating  the  inner  equations 
as  described  in  section  8.  . Since  y occurs  on  both  sides 

of  the  equation,  (61)  must  be  solved  iteratively. 

The  inner  variables  and  the  normalizations  for  (61)  are 
computed  as  follows: 

q - (5q  + | In  (l-2y ) ) (62) 


r 


where , 


and. 
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2 _ 2c  6 cyo 

r " TIW 


n - -4- 


akwoV 


n - 


2fta 

V 


Y - X* 


kV 


k - 2ita 
L 


V - 


(64) 


«a2”U0  * i In (l-2i ) ) (1  + tanh  C ) (1-2S) 


“o  "l  NL 


and  from  equation  (21a) 


86 


2 1/2 


(65) 
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Results  of  tearing  mode  analysis: 

We  carried  out  the  analysis  of  the  above  equilibrium 
for  the  following  range  of  parameters: 


B.  " 1 ' 50  55  1 * XA  “ 1 * 


1.02x10"®  e - M|°~KG  , t>  « .12  y =4irxl0~ 7 

in  weber  nr 

2 

1.38X10-22  k<?  - , m.  = 3 . 2xl0-27kg 

o^  sec^  1 


.05  m , | = .05 


and. 


N_  = lxlO1®  ions  to  N_  = lxlO18  ions 

L L 


.025  to  6 = .425 


where  6 is  the  normalized  angular  velocity  defined  by  (60). 

As  we  discussed  in  the  previous  sections,  the  effect 

of  the  resistive  layer  on  the  tearing  modes  is  contained 

entirely  in  the  criteria  of  merit,  F,  whereas  the  outer 

regions  affect  the  tearing  modes  primarily  through  the 

mismatch  Finally,  there  is  a feedback  between  outer 

and  inner  regions  that  is  particularly  sensitive  to  the 

position  of  the  outer  wall,  since  the  position  of  the 

outer  wall  directly  affects  A ' and  indirectly  affects  F 
-2 

• no  <4 6 tiiSitwi-witcu  uuw  vc  * a uc k/cuuo  lu  xcquxum 

— 2 

order,  on  the  parameters  q and  A , since  for  the  equilibria 


» y 
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being  considered,  the  parameter  r < 10  while  — ^ — , accord- 
ing to  Table  1,  varies  from  0 to-. 23.  This  range  of  values 
r D 

of  — ^ — has  been  shown  to  have  only  a minor  effect  on  the 
results. 

More  important,  as  shown  in  section  8,  it  is  basic  to 
the  tearing  mode  analysis  that  the  exponent,  h,  defined  by 
equation  (21a),  be  real.  Furthermore,  the  asymptotic 

• 'V/ 

analysis  requires  that  h 'v  0.  Equation  (65)  shows  that  real 
h is  guaranteed  for  6 < .5,  while  Table  1 inicates  that 
reasonable  results  can  be  expected  from  the  asymptotic 
analysis  for  6 < .2.  However,  we  carry  out  the  analysis  for 
values  of  6 as  high  as  .425  and  expect  the  results  to  have 
qualitative  validity. 

Exact  values  of  the  growth  rate  could  be  carried  out 
using  the  following  procedure: 


j 
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1.  Fix  the  position  of  the  outer  wall  at  some  point  at 
which  the  equilibrium  has  an  unstable  tearing  mode  for  the 
above  values  of  the  equilibrium  parameters  and  angular  velocity 
corresponding  to  6 * .025.  As  mentioned  above,  the  criteria 
for  such  an  instability  is  that  A*  >0  where  A'  is  computed 
from  equation  (30b). 

2.  At  this  low  value  of  the  angular  velocity,  the  param- 
-2 

eter  A governing  the  inner  equations  is  effectively  zero. 

Thus  the  criterion  of  merit,  F,  that  characterizes  the  effect 
of  Hall  current  and  rotation,  takes  a value  approximately 
equal  to  its  vertex  value  of  F = 2.12  as  shown  in  Figure  3. 

The  normalized  growth  rate,  y,  can  then  be  calculated  from: 

- -3/5  (4')4/S 

r - n [p  J 

where  n is  computed  from  (64). 

3.  The  angular  velocity  specified  by  5 is  increased  by 
a small  amount  keeping  the  particle  number,  N , fixed.  The 

li 

outer  equations  are  reintegrated  to  obtain  a new  value  of  A’. 

The  value  of  q is  determined  from  equation  (62)  and  since  F 
_ 2 

depends  also  on  A which  by  equation  is  a function  of  the 
growth  rate  y,  the  growth  rate  must  be  determined  this  time 
from  equation  (61)  by  iteration.  This  iteration  is  initiated 
by  using  the  value  obtained  for  y from  step  2 as  the  first, 
guess. 


-49- 


4.  Step  3 is  repeated  for  the  case  of  no  Hall  current, 

i.e.,  e ■ q ■ 0.  The  values  of  the  growth  rates,  are  used  as 

normalizing  factors;  i.e.,  Y normalized  * ¥■ 

Ye  = 0 

In  this  way  a table  of  normalized  growth  rates  as  a 
function  of  angular  velocity  for  a constant  wall  position 
would  be  obtained  for  one  family  of  equilibria  parametized 
by  N_ . However  there  is  one  problem  which  makes  the  im- 
plementation  of  the  above  procedure  difficult.  The  deter- 
mination of  A*  requires  the  evaluation  of  (30b)  at  the  edge  of 
the  boundary  layer.  However,  as  shown  in  Table  1,  the  thick- 
ness of  the  boundary  layer,  ^ , as  calculated  by  equation 

(23),  varies  directly  with  angular  velocity  taking  on  values 

on  the  order  from  0 to  10-3.  However,  this  close  to  the 

neutral  point,  where  B (r  ) = 0,  the  solutions  of  the  outer 

z o 

equations  are  asymptotically  large  and  the  integration  routines 
are  no  longer  sufficiently  accurate.  It  turns  out  that  for 
some  parameter  ranges  there  exists  a small  annulus  about  the 
neutral  point  in  which  the  plasma  position  is  close  enough 
to  the  neutral  point  for  a computation  of  i'  to  be  made 
yet  for  enough  away  from  the  neutral  point  to  enable  the 
outer  equations  to  be  accurately  integrated,  but  this  is 
generally  not' the  case. 

We  can  avoid  this  difficulty  cf  integrating  the  outer 
equations  if  we  are  willing  to  settle  for  overestimates 
of  the  growth  rates  and  incomplete  information  as  to  the 


exact  position  of  the  outer  wall.  We  merely  fix  the  value 
of  A ' at  its  maximum  value  over  the  range  of  values  on  the 
interval  of  angular  velocities  of  interest.  This  corre- 
sponds to  some  position  of  the  outer  wall.  We  see  from 

equation  (61)  that  this  value  of  A'  yields  an  overestimate 

— —2 
of  y whereas  from  equation  (63)  for  A , overestimates  of 

Y yield  smaller  values  of  A for  the  same  value  of  q thus 

from  Figure  3 this  results  in  lower  values  of  F and  once 

again  from  equation  (61),  overestimates  of  y*  We  have 

implemented  this  procedure  and  presented  the  results  in 

Figures  4-6,  where  ^ is  plotted  against  the  normalized 

Ye  •=  0 

angular  velocity,  5,  for  different  but  unknown,  outer  wall 
positions  characterized  by  and  varying  directly  with  A* 
and  for  three  values  of  the  number  of  particles, 

Nl  ■ 10^®,  10^,  10^®.  Actual  values  of  the  rigid  body 
angular  velocity,  fi,  and  ion  temperature  T are  indicated 
in  figure  4. 

Figures  4-6  validate  the  previously  stated  effect, 
namely,  that  the  Hall  current  results  in  increased  stability 
to  tearing  modes  when  the  plasma  column  rotates  in  one  direc- 
tion and  decreased  stability  when  the  column  rotates  in  the 
opposite  direction.  Figures  4-6  show  that,  in  general,  both 
stability  and  instability  are  enhanced  by  rotating  the  column 
faster  in  either  direction,  although  above  a critical  value 

of  the  wall  position,  a',  the  normalized  growth  rate,  , 

Y,  = n 
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attains  a maximum  with  respect  to  the  normalized  angular 
velocity,  5,  for  rotations  in  one  direction  and  a minimum 
for  rotations  in  the  opposite  direction. 

Figures  4-6  also  show  that  the  growth  rates  are  very 
sensitive  to  the  wall  position,  a',  and  that  the  plasma 
column  becomes  more  stable  to  tearing  modes  for  outer  wall 
positions  nearer  to  the  plasma  vacuum  interface,  i.e.,  for 
smaller  values  of  A*.  Furthermore,  for  sufficiently  small 
values  of  a',  at  a critical  value  of  the  normalized  angular 
velocity,  6,  the  tearing  mode  instability  is  quenched,  i.e., 


Figures  4-6  also  show  that  for  a given  wall  position, 

A , the  stabilizing  effect  is  more  pronounced  at  lower  values 

of  the  number  of  particles,  N , and  higher  ion  temperatures,  T. 

Li 
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Conclusions ; 

1.  Overestimates  of  the  effect  of  Hall  current  and 
rotation  on  the  growth  rate  of  tearing  modes  can  be  determined 
from  characteristics  of  the  equilibrium  at  the  neutral  point 
of  the  reversed  magnetic  field  provided  that  the  equilibrium 
is  known  to  be  unstable  to  tearing  modes. 

2.  The  effect  of  the  tearing  mode  is  characterized  by 

the  criteria  of  merit  F which  depends,  to  leading  order,  on 

-2 

only  two  parameters  q and  A . When  q > 0 the  Hall  current 

always  enhances  the  stability  of  the  plasma  system  to  tearing 

modes  while  values  of  q < 0 always  result  in  reduced  stability. 
_ 2 

The  parameter  A is  governed  primarily  by  the  position  of  the 
outer  wall  and  the  value  of  the  angular  velocity. 

3.  For  a given  number  of  particles,  the  stabilizing 
effect  of  the  Hall  current  is  enhanced  by  wall  positions  closer 
to  the  plasma  vacuum  interface.  For  wall  positions  closer  to 
the  plasma  than  a critical  value,  the  tearing  mode  instability 
is  quenched  at  a critical  value  of  the  angular  velocity  of  the 
column.  For  wall  positions  larger  than  this  critical  value, 
the  growth  rate  attains  a minimum  at  one  value  of  the  angular 
velocity  for  rotation  in  one  direction,  and  attains  a maximum 
for  another  value  of  angular  velocity  for  rotation  in  the 
opposite  direction. 

4.  For  a given  wall  position,  the  stabilizing  effect 
of  the  Hall  current  is  more  pronounced  for  fewer  numbers  of 
particles,  N. , and  higher  ion  temperatures,  T. 

AJ 
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Table  1 

Validity  of  the  Analysis. 
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■ Figure  Captions 

1:  Configuration  of  the  plasma  system. 

2a:  Analytical  solution,  £(Z),  the  resistive  layer 

equations  for  a special  case  (numerical  solution 
coincides  with  the  analytic  solution  within  the 
accuracy  of  the  graph) . 

2b:  Analytic  solution,  y(Z),  of  the  inner  equations  for 

a special  case  (numerical  solution  coincides  with 

analytic  solution  within  the  accuracy  of  the  graph) . 

3:  Dependence  of  the  criteria  of  merit,  F,  on  the 

_2 

parameters  A and  q governing  the  resistive  layer. 

4:  Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  ^ on  the  normalized  angular 

Ye*0 

velocity  of  the  plasma  <S , and  a measure  of  the  outer 
wall  position,  a’,  for  particle  number,  = 1016. 

The  actual  angular  velocity  of  the  plasma  column  is 
noted  on  the  abscissa. 

5:  Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  on  the  normalized  angular 

Ye»0 

velocity  of  the  plasma  <5,  and  a measure  of  the  outer 
wall  position,  A ’ , for  particle  number,  » 1017. 

The  actual  angular  velocity  of  the  plasma  column  is 


noted  on  the  abscissa. 


Figure  6:  Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  ^ on  the  normalized  angular 

' e*0 

velocity  of  the  plasma  6 , and  a measure  of  the  outer 

wall  position,  a',  for  particle  number,  N.  = 1018. 

L 

The  actual  angular  velocity  of  the  plasma  column  is 
noted  on  the  abscissa. 
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